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We have investigated effects of the QCD phase transition on the relic GW spectrum applying 
several equations of state for the strongly interacting matter: Besides the bag model, which describes 
a first order transition, we use recent data from lattice calculations featuring a crossover. Finally, we 
include a short period of inflation during the transition which allows for a first order phase transition 
at finite baryon density. Our results show that the QCD transition imprints a step into the spectrum 
of GWs. Within the first two scenarios, entropy conservation leads to a step-size determined by the 
relativistic degrees of freedom before and after the transition. The inflation of the third scenario 
much stronger attenuates the high-frequency modes: An inflationary model being consistent with 
observation entails suppression of the spectral energy density by a factor of < 10~ 12 . 



Theories of inflation predict the existence of a back- 
ground spectrum of gravitational waves (GWs) which 
has been created in the very early universe together with 
scalar perturbations of the energy density. Since then 
the GWs propagate freely through spacetime still retain- 
ing information about the conditions under which they 
originated. This is because their cross section is very 
small: Particles with larger cross section, such as neutri- 
nos or photons, decouple much later and therefore carry 
no information about the universe at such early times 
and high energies. This advantage of GWs renders their 
direct detection impossible until today. However, ob- 
servations of the binary pulsar B1913+16 (Hulse- Taylor 
pulsar) provide strong evidence for energy loss through 
emission of gravitational radiation [1] . 

In this article we demonstrate that in spite of being 
decoupled, relic GWs can show imprints of subsequent 
cosmic events, e.g. phase transitions. We present the 
shape of the relic GW spectrum after different scenarios 
of the QCD phase transition. Before the transition, at 
about 200 MeV, the universe is filled with a fluid contain- 
ing relativistic quarks (up, down and strange), gluons, 
leptons and photons. After the transition, quarks are 
confined into hadrons of which the pions are the only de- 
grees of freedom which are still relativistic. From this 
reduction of degrees of freedom and the assumption of 
entropy conservation, the impact of the phase transition 
on the spectrum of GWs can be derived. Numerical cal- 
culations show that the shape of the spectrum after con- 
finement and chiral symmetry breaking does not strongly 
depend on the order of the phase transition. 

However, there are scenarios for the cosmological QCD 
phase transition where entropy conservation is consider- 
ably violated and the gross features of the spectrum are 
not fixed by the estimates mentioned above. We will 
discuss a model which includes a short period with dom- 
inating vacuum energy density causing inflation. The 
melting of this energy is associated with a large entropy 
release. Within this model, the energy density in high- 
frequency GWs is much more diluted than within the 
scenarios where entropy is approximately conserved. 

We emphasize that the relic GWs, generated long time 



before the QCD transition, are not the only ones which 
are of interest in our context: During a first order trans- 
ition, bubbles of the new phase form within the old 
phase and grow under certain conditions until the fluid 
is completely converted into the new phase. This pro- 
cess involves bubble collisions and turbulences which are 
a source of gravitational radiation. For a standard first 
order transition the resulting GW spectra are already 
explored (see, for example, [2-4] and references therein), 
but a careful investigation within an inflationary QCD 
transition is still to be done. 

The structure of this article is as follows: After having 
collected a few important equations in the next two sec- 
tions, we describe a model which allows for a first order 
QCD phase transition being consistent with both obser- 
vation and lattice data (Section III). In Section IV some 
analytic estimates and the results of the numerical calcu- 
lations are presented. The calculations have been done 
within the bag model, using various parameterizations 
of lattice data, and within the inflationary scenario. We 
conclude in Section V. 



I. GRAVITATIONAL WAVES WITHIN A FRW 
BACKGROUND 

In the very early universe, a period of inflation is as- 
sumed to have taken place. This leads to perturbations 
within a homogeneous background. The tensorial part of 
these perturbations consists of GWs which are described 
by the metric 

ds 2 = a 2 (dr/ 2 - (<% - hij) Ax l dx j ) . 

Einstein's equation for this metric and an energy- 
momentum tensor of an ideal fluid gives the equation of 
motion for a GW within a Friedmann-Robertson- Walker 
background: 

+ mb! l2 - Ah i:j = 0. (1) 

H = a! I a is the conformal Hubble parameter and primes 
denote derivatives with respect to conformal time. Now 
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the expansion 



h i:j (r],x) = J ^ 3 /2 eij(k)/i k (7?)e 



(2) 



with eij (k) being the polarization tensor of the mode k, 
is applied to (1): 

hUv) + ZUKiv) + k 2 h k ( v ) = 0. 
With Vk(r)) := ahk(rj) we get 
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(see [5] for comparison). We solve this equation in order 
to calculate the behavior of the GW spectrum during 
the QCD phase transition. Therefore we need a" /a as a 
function of time or temperature: 
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a 2 (p - 3p) . 



where the Fricdmann equations have been used. We see 
that for our calculations we just have to know the evolu- 
tion of both the scale factor a and the trace of the energy 
momentum tensor. 



II. THE BAG MODEL 

A description of confinement and asymptotic freedom 
is given by the bag model [6] which assigns a positive, 
constant energy density B to the interior of hadrons. At 
very high densities or temperatures, the hadrons might 
overlap forming a single bag in which the quarks and 
gluons can move almost freely. In this way, the bag model 
can describe the confincment-deconfinement phase trans- 
ition of strongly interacting matter. 

Assuming a gas of freely moving relativistic quarks and 
gluons above a critical temperature T c and a gas of freely 
moving relativistic hadrons below T c , one can write down 
the entropy density for both phases: 



2tt 2 

s(T) = —g(T)T 3 



with g(T) := 



if T > T c 
if T < T c . 



gi and g^ are the numbers of effective relativistic de- 
grees of freedom. For g\ = 51.25 we take into ac- 
count gluons, up and down quarks, photons, electrons, 
positrons, (anti-)muons and neutrinos; 52 = 17.25 only 
includes pions instead of quarks and gluons. The Max- 
well relation 
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Figure 1. (Color online) Sketch of a possible phase diagram for 
strongly interacting matter. The early universe is generally 
assumed to evolve at very small baryon chemical potential, 
where lattice calculations predict a cross-over. At higher 
a first order phase transition is possible. The evolution of the 
universe is indicated by arrows, see also [7]. 




if T > T c 
if T < T c 

if T > T c 
if T < T c 



(4) 
(5) 



where the bag constant is fixed by T c and the difference 
in the degrees of freedom in both phases: 
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III. 



AN INFLATIONARY QCD PHASE 
TRANSITION 



Lattice QCD calculations show a crossover between the 
hadronic phase and the quark-gluon plasma at temper- 
atures of about 150-200 MeV. However, until now LQCD 
makes reliable predictions only for small chemical poten- 
tial. For example, at chemical potentials of order /ib ~ T 
the character of the transition between the two phases 
is unclear. In principle, a first order phase transition is 
possible in this region of the QCD phase diagram. 

Within the standard model of cosmology the baryon 
chemical potential during QCD time is expected to be 
well within the cross-over region of the phase structure, 
leading to the assumption that a cosmological first or- 
der transition is unlikely (see Figure 1). More precisely, 
the baryon asymmetry calculated from CMB and BBN 
data is 77b ^ 6 • 10 -10 [ ], and for small values of r/b the 
thermodynamics of a relativistic Fermi gas gives 



yields 



5gq Mb 
77r 2 q T ' 
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where g is the effective number of helicity states and in 
g q only quarks are accounted for. There are nevertheless 
different ways to circumvent this scenario and to pre- 
serve a first order transition: It has been pointed out in 
j )] that a high lepton asymmetry could be sufficient. An- 
other possibility is considered in [10]: If a short period 
of inflation occurs during the QCD phase transition, the 
baryon density could be diluted such that the correspond- 
ing chemical potential starts at values ^(Tqcd) ~ Tqcd 
(allowing the transition to be first order) and nonetheless 
matches the observations today. For example, the high 
baryon density could be provided by an Affleck-Dine ba- 
ryogenesis before QCD time [11]. The possibility of an in- 
flationary QCD phase transition has also been considered 
in [12, 13]. [14] discusses a transition with at least strong 
supercooling. 

We will be concerned with the effect of such an inflation 
on the relic GW spectrum. We only mention that also 
structure formation is influenced, primordial magnetic 
fields can be produced, and new GWs can be generated 
through bubble collisions and turbulences, [10, 15]. Thus, 
the scenario describes a fluid with high baryon asym- 
metry (see Figure 2) which is cooled down and reaches a 
quasi stable state still having a considerable amount of 
potential energy compared to the ground state. This po- 
tential energy plays the same role as the bag constant B 
in the bag model. However, if the transition into the true 
vacuum state occurs late, the universe will be dominated 
by a vacuum energy leading to inflation. Baryon dens- 
ity and thermal energy density are both strongly diluted. 
Then the vacuum energy melts and energy conservation 
gives the temperature after reheating. Afterwards, the 
fluid returns to the conventional evolution. We will im- 
plement this scenario in a simple way, not requiring a 
field theoretical description. 



IV. CALCULATIONS AND RESULTS 



A. Horizon Entry 



200- 

175- 

ISO- 

>125- 
o 

s 

g 100- 

50- 
25- 



.crossover 



fn/T-1 



Quark Gluon Plasma 
Chirally symmetric phase 



Hadron Gas 
Chirally broken phase 



color super- 
conducting 
"today" \ P hases 



200 400 600 800 1000 1200 1400 

Baryon chemical potential [MeV] 

Figure 2. (Color online) Possible track of the background fluid 
through the phase diagram for strongly interacting matter. 
Within this scenario the system intersects the first order line 
starting from high chemical potential. It is supercooled during 
inflation after which the actual phase transition takes place 
and reheating sets in. Afterwards, the evolution follows the 
commonly accepted path. 
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Superhorizon modes (krj -C 1) of relic gravitational 
radiation have a flat spectrum: 



Vt upeI (k) oc 



\h u \ 2 k 3 oc k°. 



(6) 



This holds true independently of changes in the equation 
of state. However, if the universe is not inflating, modes 
which have exited the horizon during inflation reenter 
the horizon again. This leads to a fc-dependent power 
spectrum because horizon entry occurs at different times 
for different modes and inside the horizon the modes start 
decaying: h oc I /a. For a radiation dominated universe, 
we have 



a oc rj 



v oc exp(±ikr]) . 



Figure 3. (Color online) Behavior of a GW entering the ho- 
rizon. For this calculation the background was assumed to 
be purely radiative. A phase transition is not included. The 
frequency of the mode was chosen to be the critical frequency: 
k = 2-kv* . The dashed line indicates a function proportional 
to 1/a fit to the amplitude of the GW after horizon entry. 



(see Figure 3). We denote the scale factor at horizon 
entry of the mode k with ai n (k) and see: 



Oin(fc) 



k 



oc kaf u (k) 



ain(fc) oc i. 

k 



H in {k) 

Therefore, the power spectrum after horizon entry is 
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p^ {k ) = P^(k)(^^j oc^. (7) 

Inserting the expansion (2) into the expression for the 
energy density of gravitational waves [~>] , 
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((dohijXdotii)), 



we get the fc-dependence of their energy fraction per log- 
arithmic frequency interval: 



n g (*) 



' dPs ock 5 \h k \ 2 ock 2 V g (k). (8) 



p crit dink 



Q g does not depend on k after horizon entry during ra- 
diation domination. 



B. A Step in the Spectrum 

The QCD phase transition imprints a step into the 
primordial differential energy spectrum of GWs Q, s {k). 
Let us first estimate some characterizing properties of 
this step. For comparison see the discussion in [16]. 

The typical frequency scale. We calculate the fre- 
quency of modes which enter the horizon at the end of 
the transition. The Hubble parameter for the critical 
temperature is 
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P(T C ) 
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92 T c 4 , 



3 r " ' 3 30' 
which corresponds to a comoving wave number k* with 

a c H c 



k* 
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In order to calculate the corresponding frequency today 
we need the ratio of the scale factors today (a ) and at 
QCD time (a c ): 
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To ( go 



1/3 



So for the physical wave number today, k$ , we obtain 
fc * = ^(9oV /3 



Then the frequency today is 



3.53 



180 McV 



I 92 \V 
V 17.25/ 



10" 9 Hz. 



For this result the values Tq = 2.35 • 10~ 4 eV and 
go = 3.91 have been inserted. In spite of being nonrc- 
lativistic today the neutrinos contribute to go because 
they do not transfer their entropy to the photons after 
decoupling. Thus, there is no shift in the standard 



redshift relation T cx a . 

The step-size. In order to show that J7 g is damped 
for high frequencies by a factor of 0.7, we collect the 
information of equations (6), (7) and (8) to obtain for 
the fc-dependence of 57 g 

Q g (fc) cx k 2 a 2 a {k). 

Next we use entropy conservation and get for radiation 
dominated periods: 



H cx T 2 , 



r- y/3 -i/6 -2 



g2/3 a 2 



and replace k by the Hubble length at horizon entry (k 
H in (k)a in (k)): 



fig(&) cx H 2 n af n cx g k 
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With gk we denote the number of relativistic degrees of 
freedom at horizon entry of the mode k. Assuming a 
radiation dominated universe well before and after the 
QCD phase transition leads to 



f2 g (V 3> v* 



1/3 



0.696. 



The slope of the step. Let us revisit the question of 
how the energy spectrum looks like after horizon entry. 
We will do the calculation for an arbitrary equation of 
state p = wp with the only restriction that w(a) varies 
slowly with a. We start from energy-momentum conser- 
vation, 



din If (a) = 
From this we get 

d In v = d ln(ffj n ai, 

which we use in 

din f^g(^) 
dlnu 

We write 

'„2„2 



-(1 + w(a))dlna. 



-(1 + 3w(a))dlna 
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(9) 



-|(l + 3w) din, 
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noting that the proportionality is only strict if uu in (9) 
does not depend on a. So 



dln£l(is) — 3w(cii n ) 



din v 



1 + 3w(a in ) 



(11) 
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C. Bag Model 

To keep things as simple as possible we first employ 
the bag model for a numerical calculation. So we assume 
a pressure and energy density according to equations (4) 
and (5), which leads to 



p{T) - 3p(T) = 



AB if T > T r 







if T < T c 



(see Figure 4). Note that within this setup, we do not 
account for supercooling: We assume a phase transition 
beginning exactly at T = T c . Correspondingly, we as- 
sume that entropy is conserved during the transition. 

The corresponding calculation has already been done 
in [16]. Our result can be seen in Figure 7. As expected 
from our estimates, the energy density of relic gravita- 
tional waves is reduced by a factor ^0.7 for modes with 
frequencies higher than the critical frequency scale v* . 



D. Lattice Data 
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Figure 4. (Color online) The trace anomaly of strongly in- 
teracting matter, normalized to the temperature. Within the 
bag model, the only contribution comes from the bag constant 
B. The dashed curves depict parameterizations of results of 
lattice calculations. More explanation is provided in the text. 



In order to build on a more realistic description of 
strongly interacting matter we take into account results 
from LQCD in this section. In particular, we use data 
published by the Bielefeld-BNL/RIKEN-Columbia col- 
laboration in [17] which are based on calculations includ- 
ing three quarks (u, d, and s) with physical strange quark 
mass and two degenerate light quark masses being one 
tenth of the strange quark mass. The equation of state 
has been calculated for two different improved staggered 
fermion actions, the asqtad and p4 actions. They are 
both 0(a 2 ) where a is the lattice spacing. The published 
results are mostly extracted from calculations on lattices 
of size 32 3 x 8 where N — 32 is the spatial extent and 
N T = 8 is the temporal extent. We have also calculated 
with data published by the Wuppertal-Budapest collab- 
oration in [18, 19] which are obtained using a different 
type of staggered fermion action, the stout action. 

The only things we use as an input in our calculations 
are the parameterizations of the trace anomaly given in 
[17, 19]. In order to better account for the low temper- 
ature physics, data from hadron resonance gas calcula- 
tions have been incorporated over a range of temperat- 
ures 100 MeV < T < 130 MeV. They slightly reduce the 
peak in the transition region. The label "HRG" is given 
to the quantities resulting from this procedure. In the 
"stout" case, the parameterization is chosen such that it 
also reflects the HRG results in this region. 

The pressure of the medium is given by 

P(T) _p(T Q ) = f T dT , P-3p 

T 4 T 4 T , 5 

which, together with the parameterizations of the trace 
anomaly, also fixes the energy density. The temperature 
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Figure 5. (Color online) The pressure of strongly interacting 
matter, normalized to the temperature. The dashed lines 
correspond to lattice data. For bag model 1 we included pions 
in the low temperature phase and we excluded strange quarks 
from the high temperature phase: g qU ark s = 37, ghadrons = 3. 
For bag model 2 it is the opposite way around: <? qU arks = 47.5, 
ffhadrons = 0. The lattice results tend to the Stefan-Boltzmann 
limit of bag model 2. 



dependence of both quantities is shown in Figures 5 and 
6. 

Knowing the energy density and the trace anomaly as 
functions of temperature we can now calculate the spec- 
trum of relic GWs after a QCD transition described by 
lattice data. The result displayed in Figure 7 looks very 
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Figure 6. (Color online) The energy density of strongly inter- 
acting matter within the bag model and according to lattice 
data. More information is provided in the caption of Figure 
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Figure 7. (Color online) The energy density of GWs per 
logarithmic frequency interval after different types of QCD 
phase transition. Modes with high frequency are damped by 
a factor of ~ 0.7 compared to low frequency modes. The 
solid line represents the result of a bag model calculation with 
T c = 180 MeV. The dashed curves are calculated using lattice 
data. In three cases also the hadron resonance gas ansatz has 
been taken into account. 



be fully relativistic. However, strange quarks and pions 
have masses of about the temperature scale at the trans- 
ition, so we certainly oversimplified the situation with 
our assumptions. 

This discussion reveals that only if we restrict ourselves 
to calculate solely the effect of the QCD phase transition 
with given g\ and 32, the step-size is given by 0.7. If 
we include the evolution before and after the transition, 
every loss of relativistic degrees of freedom heightens the 
step. For example, including relativistic strange quarks 
and letting the pions become completely non-relativistic 
leads to a step of size 



^without \ 
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61.75 
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= 0.613. 



These extreme cases are chosen in Figures 5 and 6 where 
we see that the lattice results lie in between. Note that 
in these diagrams only strongly interacting particles are 
considered. The total energy density and pressure are 
bigger by the amount of photons and leptons. 

After all we should state that the step-size calculated 
with lattice data depends on the temperature of the sys- 
tem at the time when the calculation ends (T 2 ). The 
displayed p4 and asqtad curves are obtained for T 2 ~ 
50 MeV. As T 2 is varied from 10 MeV to 100 MeV, the 
damping factor sweeps out a range of less than one per- 
cent (for the p4 and asqtad parameterizations) and of 
about two percent, respectively (for the p4+HRG and 
asqtad+HRG parameterizations). Unfortunately, the in- 
terpretation of our findings is not quite as straightforward 
as claimed above: The parameterizations for the trace 
anomaly in [ ] diverge for T — > and certainly lose reli- 
ability below 50 MeV. Maybe this is the case well above 
that temperature: Even in the HRG case, no calcula- 
tion at temperatures T < 100 MeV is taken into account. 
At 100 MeV, though, the trace anomaly is still not zero, 
which implies that strong interactions are not negligible. 
For lack of data in this low temperature region, we de- 
cided to extrapolate the given parameterizations to T — 
using a simple power law dependence. At least, this 
guarantees radiation domination far away from the QCD 
scale. 

The "stout" parameterization reflects the HRG results 
also for T < 100 MeV and yields reliable results down to 
zero temperature. In this case, we ended our calculation 
at T 2 — 10 MeV. (The step-size stays the same at least 
down to T 2 — 1 MeV; however, the computing effort for 
a whole spectrum would have been much bigger.) 



similar to the one we have found within the bag model: 
In both calculations we find a step at similar frequencies 
with almost the same step-size, which is what we have 
expected. The step-size is not exactly the same: This is, 
because in the bag model calculation we have assumed 
gi = 51.25 and <? 2 = 17.25 degrees of freedom before and 
after the transition, respectively. Doing so we have taken 
the strange quarks to be non-relativistic and the pions to 



E. Inflationary Scenario 

Since, until now, we have no field theoretical imple- 
mentation of the inflationary scenario at hand, we return 
to the bag model and introduce the important features 
by hand: We start with the same temperature depend- 
ence of energy density and pressure as in the bag model 
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calculation, 



p(T) = -g 1 T* + B, 



P(T) 
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and let the system evolve without undergoing a trans- 
ition. In this way we account for the demand that it is 
trapped in a false vacuum which is characteristic for a 
first order phase transition. Below the temperature 
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a period of inflation sets in. During this time the tem- 
perature falls proportional to I /a, further supercooling 
the system. After some e-folds, the thermal part of the 
energy density is negligibly small in comparison with the 
vacuum contribution B which we then force to decay into 
the conventional particle spectrum at about 160 MeV: 
As in the bag model calculation, we take into account 
92 = 17.25 relativistic degrees of freedom. The temper- 
ature after this reheating is given by the bag constant 
and the degrees of freedom, 



T r = 



30 

327T 2 
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It lies above the temperature when inflation started, but 
below the critical temperature T c : 



30 

92 7T 2 
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(gi -52)tt ; 
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B I = T c , 



which is important if we want to avoid a transition back 
to the deconfined phase. The last statement holds true 
if <?2 > 5i/4. This relation is fulfilled even if the pions 
are not included after the transition. In our case, the 
numerical value of T r is 162 MeV. After reheating we let 
the system evolve without further manual input. Figure 
8 illustrates the evolution of energy density and pressure 
during the subsequent periods described above. 

The result of this procedure can be seen in Figure 
9. The inflationary period dilutes the energy density in 
GWs much more than an ordinary phase transition. 

For modes with v < 10 -9 Hz or v > 10 -5 Hz, r2 g (^) is 
constant. This is because the modes with v < 10 Hz 
stay outside the horizon until the radiation dominated 
era after the little inflation, whereas the modes with 
v > 10 -5 Hz enter the horizon already before the little 
inflation and stay inside until the end of the calculation. 

The modes with 10~ 9 Hz < v < 10~ 8 Hz have 
wavelengths of order the Hubble horizon when the short 
period of inflation begins. They slightly dip into the ho- 
rizon and are the first to be driven out again where they 
experience no damping anymore. So they are much less 
suppressed than the modes with higher frequencies that 
are already well within the horizon when inflation begins 
and need more time to reach superhorizon scales. During 
our calculation, the intermediate modes with frequencies 
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Figure 8. (Color online) Energy density and pressure during 
a QCD phase transition with a short period of inflation. The 
total energy density decreases proportional to a -4 until the 
bag constant B ~ (250 MeV) 4 is no longer negligible. When 
the vacuum contribution is dominant, the total energy dens- 
ity does not depend on the scalefactor a. The same is true 
for the pressure which is negative during inflation: p — —B. 
The radiative energy density becomes negligible soon after 
inflation begins. 
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Figure 9. (Color online) The energy spectrum of GWs after 
a QCD phase transition with a short period of inflation. The 
inflationary phase leads to a much stronger dilution than an 
ordinary bag model transition. 



of roughly 5 • 10~ 9 Hz < v < 5 • 10 -6 Hz cross the horizon 
three times: We set the initial conditions when all modes 
have wavelengths A > 1/H; then the intermediate and 
the high frequency modes enter the horizon in the radiat- 
ive universe before the QCD transition; the short period 
of inflation stretches the intermediate modes again to su- 



perhorizon scales. After reheating, the universe is radi- 
ation dominated again, letting the Hubble horizon grow 
faster than the wavelengths: The calculation is contin- 
ued until all modes are inside the horizon again. In order 
to estimate the final shape of the power spectrum V g at 
intermediate frequencies, the effects of all three horizon 
crossings must be taken into account: 

radiation^ -t}(1) inflation^ n-y(2) radiation^ ny{2) 



super 



super 
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(12) 

After the thermal contributions to p and p have become 
small compared to the bag constant B, the equation of 
state of the medium is p = —p. This leads to 



a ex (k) 



k 



during the little inflation. We denote the scale parameter 
at the first horizon entry with a$ and the one at the 
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Figure 10. (Color online) Temporal evolution of two modes 
during a QCD phase transition with a period of inflation. The 
mode with frequency v\jv2 is at maximal/small displacement 
when leaving the horizon. Therefore, the former appears as a 
maximum in the spectrum and the latter as a minimum. 
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With ainfi we denote any scale parameter for which the 
considered mode is still inside the horizon after the first 
horizon entry. 

Since f2 g (fc) := k 2 V g (k), we obtain 



din fl g (v) 
rfln^ 



= -4, 



which is in agreement with our numerical result plotted in 
Figure 9. In our calculation inflation lasts until the scale 
factor has grown by a factor of 1000. Since the Hubble 
parameter H is constant during that time and the phys- 
ical wavelength of the modes is stretched proportional to 
a, the modes being driven out the horizon have a fre- 
quency range of three orders of magnitude. This leads to 
an overall damping of the high frequency modes in f2 g by 
a factor of (10 3 ) -4 = 10~ 12 , affirming again the numer- 
ical outcome. This result can also be obtained by noting 
that 



fig(a) oc a 



(14) 



because gravitational radiation is rcdshiftcd in the same 
way as light. 

Now we should discuss the origin of the oscillations in 
fig for the modes with intermediate frequencies: When a 



mode exits the horizon, it stops oscillating. If this occurs 
while hk(t) w 0, the mode is frozen at a very small value. 
If this occurs while hk(t) has a maximum, the mode is 
fixed at a large value. This phase effect is similar to the 
one discussed in [20] for scalar perturbations. We illus- 
trate the effect of an inflation on GWs in Figure 10. We 
expect a related mechanism to cause the slight oscillation 
in the spectrum after a normal bag model phase trans- 
ition (see Figure 7). 



F. Towards a Field Theoretical Implementation 

We would like to extend our studies of an inflationary 
scenario for the QCD phase transition: In order to ap- 
proach our aim of a field theoretical implementation, we 
introduce a scalar field \ whose potential energy replaces 
the bag constant. During the transition, the field passes 
from a false vacuum state into the true vacuum with zero 
potential energy. This corresponds to setting the vacuum 
energy density B after the transition to zero. 

We calculate within a strongly simplified model: Be- 
fore the transition, we take the dilaton field to be con- 
stant at a value associated with high potential energy, 
i.e. we calculate within the bag model. After the short 
period of inflation, we let the potential energy transform 
into kinetic energy (kination). Doing so, we assume a 
transition to a state where the dilaton oscillates within a 
potential well with negligible potential energy. Since 



1 



P : 



V(x) 



and 



1 



P : 



x 2 



v(x), 



we get an equation of state with w = 1 =>■ poca 6 . 
However, most probably the oscillating dilaton field is 
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suppressed even much faster because, during preheating, 
it decays into other fields on timescales which are very 
small compared to the Hubble time during the QCD era 
(~ ICC 5 seconds). In this case, the "oscillation domin- 
ated" period would last too short a time in order to be 
discernible in the GW spectrum. In our calculation, wc 
assume that the oscillation of the dilaton is only damped 
by Hubble expansion over a period of time which is long 
enough to be resolved in the spectrum. Afterwards, the 
remaining energy is converted into a suitable spectrum of 
relativistic particles, in the same way as it has been done 
within the previous calculation. This scenario allows us 
to get a better feeling for the signature of a sequence of 
periods dominated by different kinds of energy densit- 
ies. Since, once more, reheating is put in by hand, we 
can freely decide when it should occur. The measure will 
be the ratio by which the energy density of oscillation is 
diluted. The sequence (12) is now modified as follows: 

q-y(l) radiation^ n-y(l) inflation^ q-y(2) oscillation ^ q~)(2) 

' super ^ ' sub ^ ' super ^ ' sub" 

or radiation 

From equation (10) we have 

<c K 7% (15) 

during oscillation domination. This leads to a power 
spectrum 

^subC^' °) = draper W ^2 £5 

as can be seen from equation (13). We therefore expect 
the energy spectrum to be proportional to i/~ 3 . If the 
oscillation dominated period lasts until the modes with 
frequencies v < ICC 9 Hz enter the horizon, the shape of 
the low frequency spectrum follows from equation (11): 
fig(^) oc v. These expectations are confirmed by our nu- 
merical results shown in Figures 11 and 12. Figure 11 
displays the spectrum of primordial GWs after an infla- 
tionary QCD phase transition with an oscillation domin- 
ated period. It lasts until the energy density is diluted 
by a factor ICC 4 = 10~ s , where, for later use, we have 
defined the parameter s to be the negative decadic logar- 
ithm of this factor. Correspondingly, the scale parameter 
grows by a factor 10 s / 6 = 10 2 / 3 . As can be seen from 
equation (15), the modes which enter the horizon dur- 
ing this time cover an interval corresponding to a factor 
10 s / 3 = IO 4 / 3 . In Fi gure 11, this interval roughly lies 
between ICC 7 Hz and 3 • ICC 4 Hz, where the spectrum 
shows the expected behavior. The second horizon entry 
of the lower frequencies happens during radiation dom- 
ination and, consequently, the resulting spectrum is pro- 
portional to v~ A . In Figure 12 we compare GW spectra 
corresponding to scenarios with various durations of the 
oscillation dominated period. The result without oscil- 
latory period is also displayed. The five further curves 
result from computations including oscillation dominated 
periods with durations given by the respective value of s. 
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Figure 11. (Color online) The energy spectrum of GWs after 
a QCD phase transition including a little inflation followed 
by a period which is dominated by the kinetic energy dens- 
ity of a scalar field. The spectrum of modes which reenter 
the horizon during this stage falls like v~ 3 , whereas the en- 
ergy density in modes which reenter during the subsequent 
radiation dominated epoch is proportional to is~ 4 . 
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Figure 12. (Color online) Energy spectra of relic GWs after 
different types of QCD transition. The parameter s indicates 
the duration of the oscillation dominated period within the 
corresponding scenario. In the case s = there is no oscilla- 
tion domination included. 



In the cases s = 12 and s = 15 we also see the behavior 
fig(^) cx v for the modes which did not enter the hori- 
zon before the little inflation. This is consistent with the 
analytical estimate (11) which was obtained from energy- 
momentum conservation and the fact that GWs decay as 
1 1 a after horizon entry. The slope of the spectrum is pos- 
itive, because the rate of entering modes is higher than 



10 



during radiation domination: This leads to a stronger 
damping of low frequency modes since they enter the ho- 
rizon earlier. 



V. CONCLUSIONS AND OUTLOOK 

We have considered the gravitational background ra- 
diation assuming a flat spectrum before the cosmological 
QCD phase transition and calculated the spectrum af- 
terwards. A few different scenarios have been taken into 
consideration: First we calculated within the bag model; 
then we made use of lattice data for zero chemical po- 
tential; finally, we explored the consequences of a QCD 
transition including a short period of inflation. 

We have seen that the loss in relativistic degrees of 
freedom imprints a step in the energy density spectrum 
of GWs. In the case of negligible entropy production dur- 
ing the transition, the step-size is determined by the ratio 
of degrees of freedom in the two phases. In contrast, the 
large entropy release within the inflationary model ac- 
companies a very different shape of the spectrum: Infla- 
tion much more attenuates the energy density of modes 
inside the horizon than it is expected in case of entropy 
conservation. 

We now work at field theoretical models which imply 
the inflationary scenario [ ]. Within such a model we 
could calculate the evolution of energy density and trace 
anomaly in order to use it in the computation of the res- 
ulting spectrum. But there is an even more important 
purpose: A consistent model within the framework of ef- 
fective field theory certainly increases the reliability of 
the scenario. For example, an important point addressed 
in [21] is whether it is possible to trap the system until 
such time as inflation has lasted long enough to dilute 
the baryon density down to the numbers required by ob- 
servation. 

We neglect the possible production of perturbations 
by amplification of vacuum fluctuations. This is because 
their spectrum is suppressed as (Aqcd/™pi) 2 with re- 
spect to a corresponding spectrum created at energies 
around the Planck scale. But an important point is the 
GW production during the phase transition by bubble 
collisions and subsequent turbulences within the plasma. 

In Figure 13, we have displayed possible GW spectra 
after a first-order QCD phase transition. The charac- 
teristic amplitude h c is defined from the energy density 
per logarithmic frequency interval (which we previously 
used) as 

0.9-10- 18 Hz h r— 
hc{v) ■■= — ^n g {u), 

where h is the Hubble parameter today in units of 
100km/(s • Mpc). The upper panel shows the character- 
istic amplitude of GWs generated during inflation and 
their spectrum after an additional inflationary period at 
QCD time. Likewise, the result of a calculation includ- 
ing both, inflation and oscillation domination is presen- 
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Figure 13. (Color online) GWs after the QCD transition and 
their detectability. The (predicted) sensitivities of three GW 
detectors are displayed. PPTA (Parkes Pulsar Timing Array) 
already constrains the GW production during the QCD phase 
transition. The corresponding spectrum is displayed on the 
basis of [3] (Kamionkowski et al.) and [1] (Caprini et al.). 
See also [22, 23]. For the amplitude of relic GWs, we use the 
largest amplitude consistent with the COBE constraint given 
in [24]. The solid line in the upper panel sketches the result 
of a calculation including a little inflation and a subsequent 
oscillation dominated period with s — 15. 

ted. Assuming an extreme value of s = 15, we might 
hope that even the spectrum of primordial GWs could 
be measured by SKA. However, a scenario with a dilaton 
field being stable for more than a Hubble time seems to 
be far from realistic. 

The GWs which are produced during the transition 
can be seen in the lower panel of Figure 13. The domin- 
ant contributions come from bubble collisions and from 
the turbulences thereafter. The spectra displayed here 
are obtained from parameterizations given in [3, 4]. They 
show the maximal GW production consistent with PPTA 
measurements. Reference [4] assumes the presence of 
magnetic fields in the turbulent plasma. The most im- 
portant deviation between the spectra is the different be- 
havior in the high frequency regime: [3] includes multi- 
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bubble collisions and therefore obtains a flatter spectrum 
a for high frequencies. This result is also confirmed 
in the more recent work [2] where, however, turbulences 
are not accounted for. Neither one of the parameteriz- 
ations includes an inflationary period. A corresponding 
spectrum still has to be calculated. 

We also display the sensitivity of three relevant GW 
detectors [22]: PPTA (Parkes Pulsar Timing Array) is 
already running and has constrained the spectrum eman- 
ating from the transition. SKA (Square Kilometre Array) 
and LISA (Laser Interferometer Space Antenna) will fur- 
ther push down the limits for GW production at QCD 
time. In the more remote future, even some features 
of the shape of the relic spectrum could be revealed by 
combining measurements of Planck [25] at v ~ 10~ 18 Hz 
and of the proposed Big Bang Observer (BBO [26]) at 
v ~ 1 Hz. Given that the Planck mission succeeds in 
finding imprints of relic GWs, BBO could potentially dis- 



criminate between the standard QCD cross-over and an 
inflationary phase transition. 
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